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Abstract 

We derive the prediction of the standard model of superluminal radio sources for the apparent 
transverse velocity of a radio source located at redshift z. The apparent velocity of the source 
is reduced by a factor of 1 + z compared to that of a similar nearby source. The cause of this 
reduction is recession of the distant source due to the expansion of the universe. The apparent 
velocity of a source can be estimated from its redshift and proper motion using the values of the 
Hubble constant and the mean densities of different energy components in the universe. We derive 
an expression for the velocity valid for the currently favored cosmological model: a flat universe 
with a nonzero cosmological constant. 
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I. INTRODUCTION 



In 1966, Martin Rees predicted that "an object moving relativistically in suitable direc- 
tions may appear to a distant observer to have a transverse velocity much greater than the 
velocity of light." 1 A few years later, such motion was discovered in very distant astronom- 
ical radio sources such as radio galaxies and quasars. The motion and the sources where it 
took place are called superluminal, that is, faster than light. The discovery was the spec- 
tacular result of a new technique known as very long baseline interferometry (VLBI). This 
technique has enabled the mapping of the morphologies of radio sources to accuracy better 
than milliseconds of arc. 

Many radio galaxies and quasars contain in their nuclei compact sources of radio emission 
with several components that appear to move apart in successive VLBI images. Their 
apparent transverse velocity of separation often exceeds the speed of light. Superluminal 
motion has been observed in over a hundred of sources. 2 However, it is no unique to radio 
galaxies and quasars. There also is a class of recently discovered galactic superluminal 
sources called microquasars. All of these sources (galactic and extragalactic) are thought to 
contain a black hole, which is responsible for the ejection of mass at high velocities. 

Superluminal motion does not contradict special relativity. In the generally accepted 
standard model, it can be explained as a light travel-time effect. Superluminal radio sources 
can be modeled by one or more radiating "blobs," moving at a relativistic velocity away 
from a stationary "core." Imagine a blob of matter starting at the core and moving toward 
an observer very fast and nearly head-on. When the blob is at the core, it emits some light 
toward the observer. After it has moved toward the observer (and slightly to the side), it 
again emits light toward the observer. Because it is closer to the observer, this light will take 
a shorter time to travel to the observer. If we ignore this fact, then we will underestimate 
the true time interval, and so we will overestimate the speed. 

More quantitatively, the observed transverse velocity of the separation of the blob from 
the core, v a , is related to the true velocity, v, and the angle to the line of sight, 9, by 

vsm9 
1 — {v/c) cos 9 

where c is the velocity of light. The angle 9 = corresponds to motion directly toward 
the observer; 9 = n corresponds to motion directly away from the observer. Equation (1) 



2 



shows that if the motion is relativistic, v ^ c, and almost in the direction of the observer, 
separation at superluminal apparent velocities will be observed. 

The derivation of Eq. (1) can be found in many textbooks on astrophysics (see, for 
example, Refs. 3,4). To our knowledge, however, none of them stresses its fundamental 
limitation: it is valid only for nearby sources. However as we will explain, most sources 
are actually very distant. The apparent transverse velocity of a distant radio source also 
depends on its redshift, defined as the relative difference between the observed and the rest 
wavelength of the emitted light, z = (A G b s — A cm )/A cm . 5 Specifically, it is reduced by a factor 
of 1+z. 6 However, this effect is only briefly discussed, 7 and its physical interpretation appears 
to be missing in the literature. In this paper, we derive the effect from first principles and 
elucidate its physical meaning. 

The proper motion of a celestial object is the change in time of its angular position on 
the sky. To estimate the apparent velocity of an object with proper motion //, located at 
the redshift z, we need its angular diameter distance, Da- Da depends on the redshift 
as well as on various cosmological parameters. To date, all analyses of superluminal radio 
sources at large distances employ the classical Mattig formula for Da, which is valid only 
for a vanishing cosmological constant, A. 9 However, cosmological observations suggest that 
A ~ 0.7 in units of the critical density. 10 We will derive an expression in integral form for 
D A for corresponding to a flat universe. Following Pen, 11 we also obtain its analytic 

approximation. 

The paper is organized as follows. In Sec. II we derive the apparent transverse velocity 
of a distant radio source in the non-relativistic approximation. In Sec. Ill we present the 
corresponding derivation within the framework of special relativity. Finally, in Sec. IV we 
present a derivation based on general relativity. In Sec. V we discuss how to estimate the 
apparent velocity of a source from its proper motion and redshift and derive a relation for 
the velocity valid for the currently favored model of the universe. We conclude in Sec. VI. 

II. NON-RELATIVISTIC APPROXIMATION 

Equation (1) is valid for a source that is at rest with respect to an observer. However, 
extragalactic superluminal sources participate in the global expansion of the universe. Con- 
sequently, they have Hubble velocities of recession proportional to their distance. For distant 
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sources these velocities are a considerable fraction of the speed of light. In this section we 
will investigate the effect the recession velocity of a source has on its apparent transverse 
velocity, in the non-relativistic limit. 

We assume that at time t^ = a blob of matter is ejected from the core of a radio 
source at an angle 9 to the line of sight of an observer located at a distance r from the 
source. The observer notices this event at time = r/c. At time t^ = At, the blob 
has moved to a distance vAt away from the source, where v is its velocity relative to the 
source. The transverse displacement of the blob from the core is then Ay = v sin 9 At. 

(2) 

If the source were stationary, this motion would be noticed by the observer at time to = 

At+r/c—(v cos 9 At) / c, because at time t^ 2 > the blob is closer to the observer by the distance 

v cos 9 At. However, due to the expansion of the universe, the source is receding from the 

observer with the Hubble velocity i>h- We first assume that both v and i>h are much less 

than the velocity of light. Then, 

/,n . r + Ar v cos 9 At , , 

f P) = At + , (2) 

c c 

where Ar = v^-At. In Eq. (2), we have implicitly used the Galilean transformation of 
velocities, that is, the blob velocity relative to the observer is a sum of the source velocity 
relative to the observer and the blob velocity relative to the source. Hence, the relation 
between the observed and true time intervals, respectively At Q and At, from the moment of 
emanation of the blob from the core, is 

At = - #) = (l + fa -(3cos9)At, (3) 

where fa = i>h/c, and the apparent transverse velocity of the blob measured by the observer 
is 

Ay v sin 9 
v = — = (A) 

a At Q 1 + (3 H - f3cos9' 1 ; 

For fin = 0, Eq. (4) reduces to Eq. (1), which is valid for a blob ejected from a stationary 

source. 

From Eq. (4) it is clear that the recession of the source (away from the observer) opposes 
the effect of the blob motion (toward the observer), in a sense that it lowers the value of 
the apparent transverse velocity of the blob. For nonrelativistic recession we have i>h = cz, 
where z is the redshift of the source. This equation yields fa = z, and hence 



We reiterate that Eq. (5) is valid only when both z and (3 are much less than unity. 



III. SPECIAL-RELATIVISTIC APPROACH 

To appear superluminal, the velocity of a blob must be relativistic. In this section we 
apply special relativity to derive a formula for the apparent velocity of a blob, moving 
away from a core located at a redshift z. According to special relativity, the Cartesian 
components of velocities measured in the reference frames O and O', moving with constant 
relative velocity V along the x-axis, are related by 

< + V IR , 

= T^vp (6a) 

where 7 = [1 — (y/c) 2 ]^ 1 ^ 2 . Unprimed (primed) quantities refer to the velocity components 
measured in the frame O (O'). We align the x-axis along the observer's line of sight to the 
radio source. The unprimed frame is the observer's reference frame; the primed one is the 
source's frame. Then, v' x = — v' cos6 and v' y = v'sin9, where v' is the blob velocity relative 
to the source. The source velocity relative to the observer is V = i>h- Hence, the Cartesian 
components of the blob velocity relative to the observer are 

t>H — l/cOS# 

Vx = l-/?'/fecos0 (7a) 
_ v' sin 6 

Vy 7 h(1-/3'/3hcos^)' [m) 

where 7H = (1-^)- 1 / 2 . 

In Eqs. (7), the time dilation effect is naturally accounted for. However, in the observer's 
inertial frame we have to take into account the extra time dilation factor that occurs because 
the distance to the emitting blob (and thus the distance light has to propagate to reach the 
observer) is changing. In the time At the emitter moves a distance v x At away from the 
observer (v x and the distance may be negative). The total observed time, At Q , is At plus 
the extra factor describing how long it takes light to traverse this extra distance (v x At/c), 

Af o = Af + ^ = (l + x )Af. (8) 

c 

Note that Eq. (8) is similar to Eq. (3), but instead of the non-relativistic approximation 
for v x , v x = t> H — v' cos9, we have used the special relativity formula (7a). The transverse 



distance covered by the blob in the time At is Ay = v y At. Hence, the apparent transverse 
velocity is v a = Ay/At = v y /(l + (3 X ). From Eq. (7) we have 

v'sinO / fa - p cos 9 \ -i 



V a = 



7h(1 -(3'(3 R cose) \ l-p'p H cos6 



+ 1- B'BvcosO) ' { ) 



or, 



Next, we have 



la (i + p u ) {\-(3' cos ey 1 ] 

7H 1 -a-^) 1/2 -^-M 1/2 = (1+ , r , 



1 + /5h 1 + /5h \1 + /3h, 
In the last step we have used the expression for the special relativitistic Doppler effect. The 

result is 

B = (12) 

For consistency of notation, we have omitted the primes on P in Eq. (12). However, it should 
be clear that in Eq. (12), as in Eq. (5), v — Pc is the blob velocity relative to the source. 

Equation (12) shows that the recession (Hubble) velocity of a radio source reduces the 
amplitude of the apparent transverse motion of its blob by the factor of 1 + z. When both 
z and P are much less than unity, the term zP cos 9 can be neglected and Eq. (12) reduces 
to its non-relativistic limit, Eq. (5). What is the range of applicability of Eq. (12)? In its 
derivation we have not made any explicit assumptions about the values of z and P, except 
that z > and P 6 [0,1]. However, when applied to the real world, special relativity 
is guaranteed to work only locally, because, in general, it is not possible to eliminate the 
gravitational field globally by a suitable choice of a reference frame. (There are no global 
inertial frames in the universe.) Thus, special relativity applies only to a limited region 
around a radio source. Consequently, though it has not been assumed explicitly, Eq. (12) is 
valid for arbitrary P e [0, 1], but only for z 1. To find its generalization for arbitrary z, 
we need to apply general relativity. 

IV. GENERAL-RELATIVISTIC APPROACH 

The metric of a homogeneous and isotropic universe is given by the Robertson-Walker 
line element: 

c 2 ds 2 = c 2 dt 2 - a 2 (t)[d X 2 + R 2 S 2 ( X /Ro)(d9 2 + sin 2 0# 2 )]. (13) 
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Here, Rq 2 is the curvature of the universe and the function S(x) equals sin(x), x, and 
sinh(x) for a closed, flat, and open universe, respectively. The function a(t) is called a scale 
factor and relates the physical, or proper, coordinates of a galaxy, r, to its fixed or comoving 
coordinates, \'- r = a X- This function accounts for the expansion of the universe; its detailed 
time dependence is determined by the Friedman equations. 12 We normalize a so that at the 
present time, a(t ) = 1. 

Photons propagate along null geodesies, ds = 0. If we place an observer at the origin 
of the coordinate system, the geodesic of the photons emitted by the source toward the 
observer is radial. The source's comoving radial coordinate is x- From the metric (13), we 
have for the photons emitted from the core 

cdt ' Fa * ,,a\ 
i, -T7T\ = d X =X- 14 

For the photons emitted later from the blob, 



*° cdt' 



X ~ Ax, (15) 



where A% is the comoving distance the blob has covered, projected on the line of sight. Its 
relation to the proper distance, A/, is Al = a(t e )Ax = (1 + z)^Ax, where z is the source's 
redshift. The proper distance is vcos9At e . Hence 

Ax = (1 + z)vcos9At e . (16) 

The subtraction of Eq. (15) from Eq. (14) yields 

Ax = 



t (1) 



(i) 



cdt' C te cdt' f to cdt' 



a(f) a(t>) J t m a(t'Y 



(17) 



Because both At e and At are very small compared to the Hubble time, Eq. (17) simplifies 
to 

cAt cAt e 
a[t ) a[t e ) 

or 

Ax + cAt = (1 + z)cAt e . (19) 

By using Eq. (16), we obtain 

At = (l + z)Ai e (l-/3cos0). (20) 



For vanishing blob velocity, Eq. (20) describes the well-known phenomenon of cosmological 
time dilation. The transverse component of the distance covered by the blob is Ay = 
vsin9At e . The apparent transverse velocity is Ay/At Q , and hence 

in agreement with Eq. (12). Thus, in this case, the prediction of special relativity turns out 
to be valid globally, that is, for arbitrary z. Why? 

The amplitude of the apparent transverse motion of a distant radio source is reduced by a 
factor that depends only on its redshift; it does not depend on the background cosmological 
model. In particular, it does not depend on the mean densities of the different energy 
components in the universe, Qa and Vt rn . Here, Q\ denotes the cosmological constant A 
expressed in units of the critical density, and Q m is the mean density of non-relativistic matter 
in the universe, also in units of the critical density. Physically, these components cause the 
universe's acceleration and deceleration, respectively. Therefore, the lack of sensitivity of 
the reduction factor (of the apparent transverse velocity) to their densities implies that the 
factor is insensitive to the acceleration or deceleration of the cosmological expansion, so its 
origin is kinematic. Mathematically, the reduction is the same for any and Q m , so in 
particular for Q\ = = 0, that is, for an empty universe. This particular cosmological 
model is called the Milne model, or kinematic cosmology. The dynamics of an empty universe 
can be completely described by special relativity, 13 which is why its prediction turns out to 
be valid globally here. However, this conclusion holds only a posteriori, that is after applying 
general relativity and finding out that the reduction factor depends only on the redshift of 
the source. 

Quantitatively, the amplitude of the apparent transverse motion of a radio source located 
at a redshift z is reduced by a factor 1 + z. The general relativistic explanation of this fact 
is the cosmological time dilation between the observer's frame and the frame of the object. 6 
The special relativistic explanation is the recession velocity of the object, resulting in the 
same amount of time dilatation. These explanations are mutually consistent, because the 
origin of cosmological time dilation is the expansion of the universe, which causes distant 
galaxies to recede from the Milky Way. We have argued that the reduction of the apparent 
velocity is in essence a kinematic effect and can be qualitatively explained within a non- 
relativistic framework, as demonstrated in Sec. II. Relativistic corrections are necessary to 
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describe the effect quantitatively. 

V. ESTIMATING THE APPARENT VELOCITY 

The physical size of an object at the redshift z that subtends the angle A<p on the sky, 
Ay, can be readily derived from the metric in Eq. (13). The result is 

Ay = D A A<p. (22) 

Hence, the apparent transverse velocity of a radio source is 

v a = nD A , (23) 

where /i = A(p/At is its observed proper motion. By measuring the redshift and the proper 
motion of a radio source and knowing the cosmological parameters, we can estimate its ap- 
parent velocity. This estimate can be subsequently used to constrain the combination of the 
parameters characterizing the source, given by the right-hand side of Eq. (21). To constrain 
only the internal parameters of the source that describe its kinematics and geometry (that 
is, v and 9), it is necessary to eliminate the dependence of v a on the redshift. Therefore, 
instead of v a itself, it is common to estimate the quantity 

v m = (l + z)v a . (24) 

Unfortunately, v m is widely called the "apparent velocity" (see, for example, Refs. 6, 
16,17,18,19). This term is misleading, because it disguises the fact that the true apparent 
velocity, v a , is affected by the recession velocity of the source. We will therefore distinguish 
between v a and v m , and will call the latter the "velocity measure." The latter is the velocity 
that an observer would measure if he/she were located at the redshift z in the vicinity of 
the source. This velocity is not the velocity we measure on Earth. We have 

v m = (l + z)fiD A = fiD, (25) 

where D — (1 + z)D A is the distance measure. We recall that an object of luminosity L has 
flux / = L/(4nD 2 L ), where 

Dl = (1 + z)D = (1 + z) 2 D A (26) 
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is the luminosity distance. 

To date, all analyses of superluminal motion in extragalactic radio sources have assumed 
that the universe has a vanishing cosmological constant at low redshifts dominated by non- 
relativistic matter. 14 For such a universe, Mattig 15 derived an analytical expression for the 
distance measure D for arbitrary fl m . Consequently, the classical Mattig formula for D 
has been used in all analyses of superluminal radio sources (see, for example, Refs. 6,16, 
17,18,20). However, observations in cosmology consistently imply that A is about 0.7 in 
units of the critical density (that is, f2 A ~ 0.7). The distance measure D (or the luminosity 
distance D L ) is sensitive to the presence of A. Thanks to this sensitivity, the Hubble diagram 
for supernovae la has been successfully used to show that the universe has a non-zero A. 
Therefore, we need a relation for the distance that accounts for the presence of A. 

Observations also strongly suggest that the universe must be very close to spatially flat 
(see, for example, Ref. 21). For a flat universe with a cosmological constant, the luminosity 
distance is (see, for example, Ref. 22) 

D L {z)=cH-\l + z) [l-n m + n m (l + x) 3 ]~ 1/2 dx. (27) 

Jo 

Here, H$ denotes the Hubble constant, and spatial flatness requires Q\ = 1 — Q m . For small 
redshifts, z <C 1, Eq. (27) reduces to D L = cH^z. The integral in Eq. (27) cannot be 
performed analytically; it is a Weierstrass elliptic function. 23 

Pen 11 derived an approximate analytic expression for Dl(z), accurate to better than 0.4% 
for 0.2 < Q m < 1, for any redshift. The value of VL m is currently known to be around 0.3. 
Pen's approximation for the luminosity distance is 

D L = cH-\l + z) [77(0, n m ) - V (z, n m )] , (28) 

where 

r](z,n m ) = 2n m 1/2 [(l + z) 4 - 0.1540(1 + z) 3 s 

+ 0.4304(1 + zfs 2 + 0.19097 (1 + z)s 3 + 0.066941s 4 ] ~ 1/8 , (29) 

and 

H— ) ■ < 30 > 

(The original notation has been slightly modified.) It is a matter of choice whether to use 
a simple numerical integral or its fairly complex analytic approximation. If one chooses the 
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latter approach, then from Eqs. (25) and (28), the velocity measure v m , or f3 m = v m /c, is 



with i] given by Eq. (29). Note that f3 m is dimensionless, as it should be. Equation (24) 
implies that 



For completeness, we provide current estimates of the Hubble constant and Vt m . From 
a joint analysis of the SDSS and the WMAP data, Tegmark et al. 24 deduced that Hq = 
70^km • s" 1 • Mpc^ 1 , and Q m = 0.30 ± 0.04 (68% confidence intervals). 

VI. CONCLUSIONS 

The apparent transverse velocity of a radio source located at a redshift z is suppressed by 
a factor 1 + z. We have derived this result within the framework of both special and general 
relativity. The underlying cause of this suppression is the recession of the source due to the 
expansion of the universe. 

Given the values of the Hubble constant and the mean densities of different energy com- 
ponents in the universe, the apparent velocity of a source can be estimated from its redshift 
and proper motion. We have derived a relation for the velocity valid for the currently favored 
cosmo logical model, that is, a flat universe with a non-zero cosmological constant. 
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